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Summary. — The latest progress in calculating electroweak gauge boson produc- 
tion in association with QCD jets at hadron colliders is summarized. Particular 
emphasis is given to the recently completed QCD one-loop calculations of W+3jets 
and Wb final states. Furthermore recent developments in improving Monte Carlo 
event generators by means of combining tree-level matrix elements with parton 
showers is reviewed. 
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1. Introduction 

Top-physics is at the heart of the Tevatron and LHC physics programme. Since its 
first observation at the Tevatron in 1995 [1, 2] at lot of effort went into measuring top- 
quark production cross sections, its mass and quantum numbers. At the LHC top-quark 
physics can offer a unique window into potential new physics at the TeV scale, see for 
instance [3]. 

When considering the semi-leptonic decay of a produced pair of top-quarks the result- 
ing final state is l^ + i?™'**-|-jets (where up to two jets might be heavy-flavor tagged). 
The very same signature is provided by a leptonically decaying W associated by a corre- 
sponding number of QCD jets. To highlight the importance of the W -\- n— jets processes 
as the dominating background to top-pair production, the LHC production cross sections 
for tt+jets and W^-l-jets are presented in Fig. 1. 

Here only a generic set of cuts, namely Et,j > 30 GeV, Et,i > 20 GeV, Ai?j(j ;) > 0.4 
and |r7j,;| < 2.5, on the final-state leptons and jets (including those from the top decays) 
has been applied. From the left panel we can infer that the probability of producing a ti 
pair in association with one or more extra jets is quite significant, an observation that is 
confirmed by the corresponding one-loop calculations [6, 7]. These production rates have 
to be confronted with the W + n-jets cross sections, displayed in the right panel. While 
the inclusive W rate exceeds the tt cross section by orders of magnitude, when asking for 
> 3 jets the rates become of same size. However, even for W + 6-jets, background to the 
semi-leptonic ti + 2-jets process, at leading order (LO) background exceeds the signal. 
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LO cross section pp -> tt + N @ 10 TeV 
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Fig. 1. - The leading-order production cross sections for fi + n-jets (left panel) and W + n-jets 
(right panel) in pp collisions at -^s = 10 TeV. For the top signal the inclusive cross sections for 
stable top quarks and for the semi-leptonic decay channel are shown. For the latter and for the 
VF+jets background a set of generic jet and lepton cuts have been applied. All cross sections 
have been calculated using Sherpa [4] employing matrix elements from COMix [5]. 



From these simple considerations it is evident that there is a strong demand for 
having predictions for the VF + n-jets processes accurate at ncxt-to-leading order (NLO) 
in QCD, reducing inherent scale uncertainties of the theoretical predictions. Furthermore 
an accurate modelhng of this class of high jet-multiplicity processes in Monte Carlo event 
generators is of major importance for the success of the ambitious LHC top-physics menu. 



2. W+3jets at next-to-leading order 

Until recently NLO predictions have been available only for final states involving a W 
boson and up to two additional jets [8]. Significant progress in the evaluation of virtual 
matrix elements involving many external legs has enabled two independent groups to 
eventually calculate W -I- 3jcts at one-loop accuracy. 

In Refs. [9, 10] the leading-color approximation to the full result has been presented. 
In this calculation the D-dimensional generalized unitarity method as described in Ref. 
[11] is used to evaluate the loop amplitudes. The actual calculation is performed in the 
framework of the MCFM code [12]. The authors of Refs. [9, 10] proposed a prescription 
called "leading color adjustment" that allows them to provide a sensible approximation 
to the full-color NLO result. In essence they rescale the leading color one-loop result 
by a constant factor defined to be the ratio of the LO full color cross section over its 
leading-color approximation. 

In Refs. [13, 14] the first complete NLO calculation of VF -f Sjets has been presented. 
This calculation includes all partonic subprocesses and is exact in the treatment of color. 
For the one-loop matrix elements the program BlackHat [15] is used that is based on 
unitarity methods [16, 17, 18]. For the generation of the real-emission matrix elements, 
the Catani-Seymour dipole subtraction terms [19], as well as all phase-space integrations 
the Monte Carlo generator Sherpa [4, 20] is used. 

From NLO calculations we can expect a reduced dependence on the unphysical renor- 
malization and factorization scales. However, they still exhibit a scale dependence. 
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Table I. - Total cross sections in pb with scale dependence for W + n-jets at Tevatron compared 
to data from CDF [21]. Numbers taken from Ref. [14], a complete description of the calculational 
setup and the cuts used can be found therein. 



In Tab. I the theoretical prediction for the W + 1,2,3-jcts cross sections calculated 
at LO and NLO are compared to a measurement by CDF [21]. The inherent scale un- 
certainties are indeed significantly reduced for the one-loop results. The newly obtained 
W + 3jets NLO result is in perfect agreement with the data. The predicted scale uncer- 
tainty for VF^ + 3jets production at the LHC is also largely reduced at NLO. Considering 
pp collisions at 14 TeV and E^^'* > 30 GeV, Ref. [14] quotes 
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Besides a reduced scale dependence of the total W + 3jets cross section the NLO 
calculation exhibits largely narrowed uncertainty bands for differential distributions. This 
is exemplified in Fig. 2, where the transverse-momentum distribution of the third-hardest 
jet at Tevatron and LHC energies is shown. However, care has to be taken which central 
scale is actually used in this intrinsic multi-scale problem. As pointed out in Refs. 
[10, 14] a choice like the bosons transverse momentum, Ej!', can yield unphysical results 
for certain distributions, originating from large kinematic logarithms. A seemingly more 
appropriate choice is the total partonic transverse energy, Ht. 
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Fig. 2. - Transverse-momentum distribution of the third-hardest jet in VK + 3jets events at the 
Tevatron, compared to data from CDF [21] (left panel) and the LHC (right panel) at leading- 
and next-to-leading order. Figures taken from [14]. 
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Fig. 3. - Jet transverse-momentum distributions in W^ + 3jets production at a 10 TeV LHC. 



The default scale is chosen to be no 



fJ-R = ^J■F 



= v^ 
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The local scale choice 



corresponds to taking each as factor in the LO calculation at its respective kx splitting scale. 
For details on the calculational setup see [22]. 



Rcf. [22] discussed the possibility to accommodate shape differences between NLO 
and LO results by appropriately choosing scales for the strong coupling factors in the 
latter. In particular a local scale choice was investigated where each as factor is evaluated 
at a reconstructed kx splitting scale. In Fig. 3 a comparison between the default scale 



fJ-o 



Pt,w 



'w 



and the local prescription for the transverse-momentum distribution 



of the three hardest jets is shown. The local scale scheme is in much better agreement 
with the NLO shapes. This approach of local as factors is commonly used in parton 
shower Monte Carlos and in particular in calculations that combine tree-level multi- 
parton matrix elements with showers [23]. The explicit comparison for W + 3jets final 
states at NLO confirms observations made in Refs. [24, 25] for W -\- l,2jets production 
and re-afBrms the predictive power of the matrix element parton shower approach. For 
a further study along these lines see [26] . 



3. W+heavy flavors at next-to-leading order 

Concerning backgrounds to top-quark production special attention has to be given to 
W^-l-jets final states with one or two jets being b-tagged. Using massive partons in the 
theoretical calculation removes corresponding soft and collinear singularities as they are 
regulated by the finite quark mass, however at the price of the fixed-order calculation 
being more difficult. At present the NLO corrections for Wbb, with massive ti-quarks, are 
known [27]. When applying cuts that suppress contributions from the threshold region 
TTijf, « 2mf, the actual difference between the fully massive calculation and the limit 
mi, = is typically less than 10% [27]. 
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Fig. 4. - LO difTerential cross section for Wb production at Tevatron as a function of the bb 
invariant mass. Only one b-quark is found inside the fiducial volume here. Figure taken from 
[28], further details can be found therein. 



However, when sensitive to the threshold region or in case that just one heavy quark 
is tagged the massless approximation is not applicable. In the latter case the unobserved 
heavy jet must be integrated over the whole phase space thus introducing reference to 
the b-quark mass, cf. Fig. 4. One way out is to use heavy-quark parton distribution 
functions - the so-called variable flavor scheme (VFS) that has the additional advantage 
to re-sum large logarithms of the type \n{mw/'mi,) to all orders. 

In Ref. [28] a full NLO calculation of producing a W boson in association with just 
a single b-jet has been presented. This calculation consistently combines the massive 
Wbb calculation of [27] with the VFS computation of Wbj [29]. This calculation is an 
important ingredient when comparing the recent CDF measurement of the W associated 
6-jet cross section [30] with the NLO QCD calculation 
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There is obviously tension between experiment and the theoretical results from NLO 
QCD as well as Monte Carlo predictions relying on matrix-element parton-shower merg- 
ing [30] . The source of this disagreement is still under study - but might be assigned to 
the scale choice in the calculations [31]. 



4. Monte Carlo event generators 

In cases we lack a full NLO calculation (e.g. W+ > 4jets) or observables are sensitive 
to multiple-parton emission and hadronization effects, theoretical predictions rely on the 
ability of multi-purpose Monte Carlo generators such as Pythia [32], Herwig [33] or 
Sherpa [4, 34] to account for the underlying physics. Over the past decade enormous 
efforts went into improving these calculations by consistently incorporating multi-leg 
tree-level matrix elements into parton-shower simulations in the spirit of [23, 35, 36]. For 
an overview of available approaches and an extensive comparison for VF-f jets production 
at Tevatron and LHC see Ref. [38] . Essentially two major problems have to be addressed 
by each tree-level merging algorithm: 
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• How to attach a parton shower to a muhi-leg tree-level matrix-element calculation 
without spoiling the logarithmic accuracy of the underlying QCD rcsummation? 

• How to avoid potential double- or under counting of phase-space configurations 
present in the parton shower and corresponding matrix-element calculations? 

To accommodate these conditions in a generic tree-level merging algorithm 

• multi-parton matrix elements get regularized through a suitably defined jet measure 
(e.g. a critical kx- or cone-like distance); 

• appropriate starting conditions for the initial- and final-state parton shower have 
to be determined and certain (hard) shower emissions need to be vetoed. 

In particular the second item is subject to certain approximations in the various 
schemes. An important concept to overcome those approximations is a so-called truncated 
shower, first proposed in Ref. [39] . The underlying observation is that due to a mismatch 
of the jet-measure, used to slice the emission phase space, and the actual shower-evolution 
variable the radiation pattern of soft/large-angle emissions can be distorted. 

In Ref. [40] such a truncated shower was implemented for the first time. The imple- 
mentation relies on the shower algorithm based on Catani-Seymour dipole factorization 
[41] and combines it with the matrix-element generators available inside the Sherpa 
framework. The method has successfully been applied to jet production in e^e~ col- 
lisions, the DrcU-Yan process [40], prompt-photon production [42] and deep-inelastic 
scattering [43]. As of version 1.2 it constitutes the default method for combining ma- 
trix elements with parton showers in the Sherpa generator. The new merging approach 
yields a largely reduced dependence on the intrinsic merging parameters compared to the 
previous CKKW implementation in Sherpa and other merging algorithms [38]. This is 
illustrated by the systematics studies for Z°/7*-|-jets production at Tevatron presented 
in Figs. 5 and 6. The first figure presents a comparison of the jet multiplicity and 
leading-jet pr distribution while in the latter the variation of the kx differential jet rates 
di_j.o and d2-^i for three different values of the slicing measure are presented. 
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Fig. 5. - Jet multiplicity (left panel) and the leading jet pr spectrum (right panel) in inclusive 
Z /7*+jets production compared to data from CDF [37]. Figures taken from Ref. [40] 
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Fig. 6. - Differential jet rates doi (left panel) and di2 (right panel) for the CDF Run II kr- 
algorithm [44]. Displayed are the predictions for three different values of the merging cut. 
Figures taken from Ref. [40] 



5. Conclusions and Outlook 

Understanding the process of electroweak gauge boson production in association with 
QCD jets is crucial for the success of the top-physics programme both at the Tevatron and 
even more so at the LHC. In the last few years there has been enormous progress in the 
calculation of one-loop corrections to multi-parton final states. As a result the processes 
W + 3jcts and Z'^ /j* + 3jets arc meanwhile known at ncxt-to-lcading order in QCD. 
In fact Ref. [45] already reports on first steps towards the calculation of ly -I- 4jets at 
the one-loop level using the BlackHat+Sherpa package. At this conference M. Worek 
reported on the Helac-Nlo package, that has proven to be capable of doing calculations 
of this complexity as well and next-to-leading order calculations for Wbb+ < 3jets now 
seem to be feasible. 

Concerning the simulation of M^-l-jets with Monte Carlo event generators a high level 
of sophistication has been reached. The approach of combining multi-leg tree level matrix 
elements with parton showers has become a widely used standard that delivers results 
in good agreement with data from Tevatron and exact higher-order calculations. One 
important future direction will be to precisely understand how these methods can be 
generalized to allow for the inclusion of one-loop matrix elements. First proposals in this 
direction have been made and implemented already, cf. Refs. [46, 47] and P. Nason's 
contribution to these proceedings. A novel procedure how to combine ncxt-to-leading cal- 
culations of different final-state multiplicity, though not facing the problem of attaching 
parton showers, has been presented in Ref. [48]. 
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